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Analytic Variational Investigation of Euclidean SU(3) Gauge Theory 
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Analytic variational techniques for lattice gauge theories based on the Rayleigh-Ritz(RR) method were pre- 
viously developed for euclidean SU(2) gauge theories in 3 and 4 dimensions. Their extensions to SU(3) gauge 
0^ ■ theory including applications to correlation functions and mass gaps are presented here. 
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1. Introduction 

Despite impressive strides made in the Monte 
Carlo(MC) studies of lattice gauge theories, a sat- 
isfactory understanding of the wavefunctionals of 
these theories is still lacking. Motivated by the 
success of variational methods of quantum me- 
chanics in addressing such issues, a program of 
developing these methods for lattice systems was 
undertaken. At first these methods were applied 
to spin systems, both discrete and continuous, in 
two dimensions, with success. Subsequently, the 
much harder problem of applying these methods 
to 4-dimensional gauge theories was undertaken 
and successfully solved for the case of U(l) and 
SU(2) gauge theories. Details of these and other 
related works along with the relevant references 
can be obtained from [EJ . 

The basic idea is to apply the RR technique 
to the spectrum of the transfer matrix, denoted 
symbolically by T. One makes an ansatz for the 
largest eigenvalue state guided by the usual intu- 
tions of nodelessness, symmetry etc. 

The simplest such ansatz for a lattice gauge 
theory is 



rp = exp(a RetrPi 



(1) 



where Pj are the relevant plaquette variables lying 
on a "time" slice with tr being taken over some 
irr. of the group, and a a variational parameter. 
Recall that T connects configurations on neigh- 
bouring time slices. The estimate for the largest 
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eigenvalue is then given by 
A = max a (ip,Tip)/(ip,ip) 



(2) 



The denominator in (2) is the partition function 
of a 3-dimensional theory with coupling 2cv, while 
the numerator can be interpreted as the parti- 
tion function of a generalised 3-dimensional the- 
ory where the links per site are doubled and with 
anisotropic couplings a + /3/2, (3. 

Our idea was to construct the appropriate 
transfer matrices for these lower dimensional 
problems and then apply the RR method to their 
spectra, and so on. Under the simplifying as- 
sumption justified a posteriori, that the same 
variational parameters are used for all the sub- 
sequent ansatze, one arrives at the variational es- 
timate for the free energy per site of the form 



F 4 = max a [F H {f3 H ) + G(a)] 



(3) 



where F H is the free energy for a 4-dimensional 
unit hypercube with coupling 0h = 3 /2a + /3/4, 
while G(a) is expressible in terms of the norms of 
the various ansatze (details are given in 0]). The 
unit hypercube partition function as well as G(a) 
are evaluated by making use of group character 
expansions 



exp{xRetr f{g)) = d r tr r (g)b r (x) 



(4) 



where trt(g) is the character in the i-th irrep. 
whose dimension is (fj. For SU(2) b r (x) = 
l2r+i(2x) / x where the I r {x) is a modified Bessel 
function. For SU(3) their analytic form is not 
known. / corresponds to the fundamental rep. 

The main result which made the analytic inves- 
tigation of the SU(2) case possible is the explicit 
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expression for the partition function of a unit hy- 
percube in terms of a product of 16 SU(2) 6-j 
symbols [Q. This result was obtained by gener- 
alising the graphical rules of Yutsis, Vanagas and 
Levinson [[|. Sybolically expressed, this result 
reads 

g _ P = rf rl b rl ....d r24 b r 24 J] 6j's (5) 

For SU(3), however, there are additional com- 
plications due to i)the complex nature of irreps, 
ii) multiplicity problems associated with weights 
and tensor decompositions. These are being in- 
vestigated at present. However, substantial con- 
tribution to (5) comes from surfaces that only de- 
pend on di and the number of singlets N rst con- 
tained in the tensor product of 3 irreps. Even 
though the latter is affected by multiplicity prob- 
lems too, a programme was developed based on 
the elegant prescription given in Q . 

In practice, the sum (5) is evaluated by find- 
ing the contribution of all possible closed surfaces 
that can be embedded in the unit hypercube. A 
complete classification of such surfaces is given 
in (|]. For example, each surface of area A and 
topology of a sphere contributes Yl^r^r to (5), 
while each surface with area A and torus topol- 
ogy contributes Yl^r ■ There are, for example, 8 
spheres with A=6 and 3 torii with A=16 in the 
unit hypercube. Likewise, surfaces obtained by 
gluing 3 disks of area A r spanned by plaquettes 
in the irrep r, A t by plaquettes in the irrep t and 
A s in irrep s, contribute J2^rstd r dtd s b^ r bf''bf t . 
More complicated surfaces require a knowledge 
of the 6-j symbols. The sum in the above formu- 
lae(except(5)) is only over nontrivial irreps. 

The b r (x) for SU(3) were calculated by gener- 
ating their series expansions upto x 100 using an 
algorithm to find the number of times the irrep 
r occurred in (/ © /)" by a recursion relation. 
The resulting b r 's were checked against a numer- 
ical evaluation of their integral representation for 
x as large as 10, and agreement was seen to dou- 
ble precision accuracy. A weak coupling expansion 
was also developed using the symbolic manipula- 
tion package FORM and again good agreement 
was seen with the series expansion. 

At the moment we have restricted the sum(5) 
to only those surfaces for which one at most needs 



the N rs t- It is encouraging(also from the point of 
more accurate extensions of this method) that in 
the range < (3 < 2.25(this corresponds to the 
range (0, 6.75) for the coupling used in literature) 
this truncation works quite well as will be shown 
by the results. 

2. Correlation Functions And Mass Gaps 

There are three classes of plaquette-plaquette 
correlations within the unit hypercube that are 
relevant for the study of mass-gaps ;these are 
a)when the plaquettes on both the time-slice are 
in the same spatial plane and b)when they are 
in planes orthogonal to one another (this class 
will be called dh3). Case a) further splits into 
al)when the two plaquettes can be embedded in 
a cube(dh2) and a2) otherwise(dhl). It is useful 
to factor out bo(x) 24 in (5) (here refers to the 
trivial irrep) and write 

-F H = 2Alnb a (x) + ln(l + A) (6) 

where A is the sum of contributions from various 
surfaces, examples of some having been given be- 
fore, except that & r 's in those expressions should 
be replaced by B r = b r /bo- The average plaque- 
tte(normalised) is then given by 

(P> =Bf{{3 H ) + \/72A'{\ + A)- 1 (7) 

The connected correlation function between two 
given plaquettes Pi and P 2 can be obtained as 
follows: let the coupling for these two be j3\ and 
@2 with all other couplings being (3h- The free 
energy in that case is given by 

e- pH =b (P H ) 22 bo(/3i)b a (/3 2 )(l + A((3 1 ,l3 2 )) (8) 

Here A contains contributions only from sur- 
faces in which both the plaquettes Pi and P 2 
are embedded. The connected correlation func- 
tion(normalised) is given by 

(PP 2 ) = l/9( J-^-)(-F ff )| /3l=/32=te (9) 

It is clear from the above that a lot fewer surfaces 
contribute to the connected correlation lengths 
than to (P). At low /3h, the dominant contri- 
bution to the average plaquette is P/(/3y), and 
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all connected correlation functions nearly vanish. 
Only dh2 receives contribution from one out of 
the 8 spheres of area 6. With increasing (3h 
other surfaces start contributing and all the three 
classes become appreciable. 

2.1. Mass Gap estimates 

Having obtained a variational "ground" state 
1 0) * for the transfer matrix, we can extend the 
technique to obtain estimates for the lower eigen- 
values and hence the mass gaps. There are some 
important technical differences in implementing 
this strategy for lattice systems as compared to 
quantum mechanics. We summarise these points 
and present the details elsewhere: a) in the vari- 
ational treatment of quantum mechanical sys- 
tems, the ansatze for the excited states can be 
chosen quite independently of the ground state 
ansatz. But in the RR applied to lattice theo- 
ries the trial excited state functionals have to be 
of the form F|0)' where F is a polynomial func- 
tional;otherwise, the mass gap grows with volume 
in the thermodynamic limit. 

b) if the excited state 'projector' F is expressed 
as F = c r 9 r with 9 r a basis for such projectors, 
a necessary condition for the problem determin- 
ing the eigenvalues of the transfer matrix to be 
well defined is 

(0|Tfl r |0)* _ (01^10)* 

(0|T|0>* (010)* 1 ' 

Otherwise, the problem of determining the eigen- 
values of ND^ 1 where the matrices N and D are 
given by N rs = ( | T | ) f and D rs = ' (0|0) ' f , 
becomes ill-defined in the sense that while D rs 
grows as volume in the thermodynamic limit, N rs 
grows as vol 2 . 

Situations where (10) is satisfied automatically 
are i) due to symmetry reasons both (OlT^lO) 4 
and (O|0 r |O)' vanish;this happens, for example, 
in the X-Y model when 6 r are rotationally non- 
invariant, ii) when 1 0) * is an exact eigenstate of 
T. iii) when |0) 4 is of the form exp(J2 \x r B r ) where 
(/x r ) are determined variationally, in which case 
(10) are the conditions for the stationarity of Ao 
w.r.t {Hr). 

In view of the last remark we should build the 
excited state projectors only out of P's. Nev- 



ertheless, in the case of SU(2), Wilson loops of 
length 6 could also be considered as it was seen 
that in the range of /3-values of interest, (10) was 
satisfied quite well. But D rs in that case required 
the knowledge of 6-j symbols; hence for SU(3) we 
have restricted our attention to only plaquettes in 
the fundamental rep. Following standard proce- 
dure, one then constructs linear combinations of 
these that transform irreducibly under the cubic 
group. With this choice, N and D are already 
diagonal. To the accuracy considered the various 
cubic irreps were nearly degenerate in mass (which 
is not the true situation) and hence only the low- 
est (scalar)mass estimates will be shown in the 
table. 

3. Results 

In the strong and weak coupling ends the anal- 
ysis can be carried out to yield a in analytic form. 
In the intermediate regions the analysis has to be 
made numerically but without MC. 

The first four terms of the exact free energy 
in both these ends are reproduced for all gauge 
groups and dimcnsions(spacc-time);actually for 
d=4 the coefft of the fourth term in the weak 
coupling is numerically very close, but not ex- 
actly equal to, the exact result. The variational 
parameter a — (3/2 + c/3 5 for small (3 and equals 
7/9/3 for large /3. 

In the case of SU(2) the average plaquette for 
d=3 was indistinguishable from high statistics 
MC; in d=4 it agreed to within .1% everywhere 
except the crossover region where the maximum 
discrepancy was about 2%. The specific heat in 
d=4 also agreed very well. The string tension and 
mass gap agreed very well upto the crossover but 
beyond there was a noticeable deterioration. For 
further details see [jjj. In the mass gap estimates 
inclusion of operators in higher irreps. made 
practically no difference while there was a dra- 
matic decrease in the mass gap when the effec- 
tive dimensionality of the transfer matrices was 
increased. Also inclusion of larger Wilson loops 
improved the estimates. 

In d=3 our variational results for SU(3) looked 
qualitatively very much like the results for SU(2) 
in d=3; there was a noticeable peak in the specific 
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Table 1 



SU(3): Local Quantities 



0H 


a 


(P) 


(P)4 


(P)an 


0.6 


.6000 


.1138 


.1168 


.1164 


1.0 


.9950 


.2028 


.2078 


.2040 


1.2 


1.1880 


.2501 


.2566 


.2517 


1.4 


1.3720 


.3012 


.3017 


.3012 


1.6 


1.5410 


.3512 


.3592 


.3522 


1.8 


1.6960 


.4019 


.4330 


.4044 


2.0 


1.8500 


.4586 


.4821 


.4554 



heat curve which however was not as sharp as the 
one in d=4.Here we did not compare the results 
with MC data. In the weak coupling region analyt- 
ical results were obtained with the help of FORM 
and the first four terms of the free energy density 
agreed with exact results. The details of this case 
have not been given here but will be reported else- 
where. Efforts are under way to enlarge the scope 
of these calculations by considering more compli- 
cated ansatze. Inclusion of fermions is also being 
considered. 

The results for SU(3) in d=4 are summarised 
in two tables. A two- fold check with MC was 
performed. The first of these consisted in simu- 
lating the unit hypercube itselfjthe 'infinite' lat- 
tice was simulated on a 4 4 lattice. One noticed 
that when (3 H was small there was no apprecia- 
ble difference between (P) obtained by the two 
simulations. But as (3 increased there was a con- 
siderable difference between them;but when the 
(P) of the hypercube simulation at (3h was com- 
pared with that of the 4-site simulation at 0, there 
was very good agreement. This can be viewed as 
a MC realisation of the method described here. 
On the other hand, since in our SU(3) studies we 
have not been able to include all the possible sur- 
faces in (5), the MC of the hypercube served as 
a check on the accuracy of the analytical results 
represented by (P) a n hi table 1 and by dh2(an) in 
table 2. The hypercube MC was performed with 
50,000 sweeps each. The tables show that as far as 
(P) is concerned, upto (3 ~ 2.25 the RR estimate 
is very good. Comparison of the connected cor- 
relations is trickier; at low /3 these are very very 
small and MC is unable to measure them reliably. 



Table 2 



SU(3): Correlations and Masses 



0H 


dh2 


dh2(an) 


mass 


0.8 


-.00002 


.00006 




1.0 


.00005 


.00018 


5.5 


1.2 


.00046 


.00041 


4.4 


1.4 


.00076 


.00078 


4.6 


1.6 


.00132 


.00121 


3.6 


1.8 


.00198 


.00148 


4.0 


2.0 


.00237 


.00129 


4.2 



At larger (3, however, the error in (P) a n becomes 
significant due to the neglect of many surfaces in 
(5). Thus reliable comparison can be made only 
in the window (1.2-2.2) which fortunately covers 
the interesting range. But even in this range only 
dh2 could be sensibly compared. 

3.1. Future Prospects 

From the behaviour of the mass gap in both 
SU(2) and SU(3) cases it is clear that the ansatz 
(1) works very well for local quantities like (P) 
but in the region close to the crossover it is not 
so efficient in describing long distance correla- 
tions. On the other hand this was the simplest 
ansatz one could have considered and it is sur- 
prising that it worked as well as it did. It is clear 
that larger operators have to be included;also the 
weakest link in the present implementation is at 
the stage where the same type of ansatz was used 
for the 2-dim transfer matrix. But in that case 
the ansatz predicts strictly vanishing d=2 correla- 
tions (connected). Any further improvement must 
circumvent this. An effective way of introducing 
larger operators is to choose T"|0)* as the trial 
ground state. This will then involve the enumer- 
ation of surfaces embedded in larger hypercubes. 
This is under study. 
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